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Adaptation and Disturbance Rejection for 
Output Synchronization of Incrementally 
Output-feedback Passive Systems 

Hongkeun Kim and Claudio De Persis 


Abstract 

This note addresses the output synchronization problem of incrementally output-feedback passive 
nonlinear systems in the presence of exogenous disturbances. Two kinds of distributed controllers are 
proposed; one placed at the nodes and the other placed at the edges. Each of them is synthesized based 
on the adaptive control method to cope with the shortage of passivity, and on the internal model principle 
to deal with the disturbances. The proposed controllers synchronize the outputs of the nonlinear systems 
when the solution of the closed-loop system is bounded. Based on this, we present a class of systems 
for which boundedness of the solutions is guaranteed. The analysis used in this note is also applicable 
to a case where systems are coupled via links modeled by dynamical systems. Simulation results of a 
network of Van der Pol oscillators show the effectiveness of the proposed methods. 

Index Terms 

Output synchronization, incremental output-feedback passivity, strong coupling, disturbance rejec¬ 
tion, solution boundedness. 


I. Introduction 

Over the past deeades, synehronization of multiple systems intereonneeted through a graph 
has received considerable attention from the control community due to its large number of 
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applications such as rendezvous, formation eontrol, frequeney synehronization of power grids, 
oseillators synehronization, to name a few. Motivated by the passive properties of physieal 
systems sueh as eleetrieal networks and oseillators, the passivity-based approaeh has proven 
its usefulness in dealing with the synehronization problem. See, e.g., |[I1-[I8| (refer to |I91 for 
other approaehes based on QUAD or eonvergent properties). In partieular, important elasses of 
nonlinear oseillators fall into the eategory of systems with a shortage of passivity for whieh 
proving synehronization still relies on passivity-based arguments ifTOll . 

When dealing with systems that exhibit a shortage of passivity, one of key requirements to 
aehieve synehronization is the so-ealled strong eoupling eondition, meaning that the algebraie 
eonneetivity of the graph should exeeed eertain threshold value El-Jll, JTl, [|9|| (even though 
there is no disturbanee aeting on the systems). Another important requirement is that the solutions 
of the elosed-loop system are bounded Q, [[6l|, flUl, whieh may not be guaranteed in general. 
This is beeause Lyapunov funetions used to show synehronization are often eonstrueted from 
the differenees between states and even when boundedness of the Lyapunov funetion is proven, 
it is not possible to infer boundedness of the states but only of their differenees. 

Many dynamieal networks are open systems that interaet with the environment (e.g., an 
inventory system with in-/out-flow of material, a smart grid subjeet to unknown demand and 
generation). As sueh, they are often affeeted by external perturbations that ean disrupt synehro- 
nization. Internal model based eontrollers have been reeently proposed as a mean to restore 
synehronization despite the effeet of disturbanees Id, IIH, [HTI . Previously, passivity-based 
internal model eontrollers were proposed to deal with formation eontrol problem with unknown 
referenee veloeity [fT2ll . Furthermore, passivity and internal model eontrollers play a role in the 
eontrol of dynamieal networks when eonstraints and optimality eonsiderations must be also taken 
into aeeount flUl, [fT3]l . lfT4ll . 

Motivated by the work diseussed above, we study the output synehronization problem of 
inerementally output-feedbaek passive (iOFP) nonlinear systems on an undireeted graph, in the 
presenee of disturbanees generated by exogenous systems. Two different struetures of distributed 
eontrollers are handled in this note. In the first one, eaeh loeal eontroller is plaeed at the 
eorresponding node, while in the seeond one, a eontroller is assoeiated to eaeh edge. Either 
at the nodes or at the edges, eaeh eontroller is a eombination of an adaptive law to eope with 
the shortage of passivity [fTSl Chapter 6] and internal model ifT^ to deal with the disturbanees. 
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It is shown that despite the shortage of passivity and the external disturbances, the proposed 
controllers enforce the output synchronization of the iOFP systems, provided that the solution of 
the closed-loop system is bounded. Then, we show that for nonlinear systems which are input- 
to-state stable (ISS) relative to a compact set X, boundedness of the solution can be guaranteed. 
A class of open-loop systems ensuring both iOFP and ISS properties is given as well. 

The structure in which controllers are located at the nodes is common in the literature and 
it is already known that (iOFP) nonlinear systems interconnected via a static diffusive coupling 
achieve synchronization under strong coupling conditions [j2| . In this note, the difficulty descends 
from the fact that as opposed to the static diffusive coupling, the systems at the nodes are coupled 
via dynamical systems that aim at guaranteeing synchronization while rejecting the action of the 
disturbances. Meanwhile, controllers at the edges have attracted interest more recently ll^. [|9]|, 
[fm . In certain problems in, e.g., distribution networks, the inputs to the dynamical systems at 
the nodes are constrained to satisfy certain physical laws (such as Kirchhoff’s laws) and having 
controllers at the edges is more convenient since they regulate the “flow” exchanged among 
the different nodes. In addition, the analysis carried out in the case of dynamical controllers at 
the edges turns out to be useful to deal with the case in which the dynamics at the edges is 
given, a feature which arises in those synchronization problems where the graph models physical 
interconnections among the nodes [fTOll . [[T3ll . [fTTll . 

A few remarks on the proposed controllers are as follows. In contrast to |j5l, (fill, each adaptive 
law assigned to the corresponding node (or edge) is one-dimensional, leading to lesser dimension 
of controllers in general. Moreover, it does not require the symmetry (presented in, e.g., [l5]|, 
(0) of the initial values and update gains of the adaptive laws. Finally, the controllers are 
driven by relative outputs of the nonlinear systems only, while additional communication of the 
partial states of the internal models is necessary in [l6j|. Thus, pure output feedback control for 
synchronization is achieved in our case. 

Notation: We denote an undirected graph by ^ = (AC, £^,A), where the node set M = 
{1,2,...,A^} is a finite nonempty set of nodes, the edge set C TV x AA is a set of pairs 
of distinct nodes satisfying (A j) & S (j, i) G E, and the adjacency matrix A = [a^] E 
is a symmetric nonnegative matrix defined in a way that > 0 if (i,j) E E and aij = 0 
otherwise. A path connecting nodes i ^ j is a sequence of distinct nodes, {pi,P 2 , ■ ■ ■ ,Pd}, 
such that Pi = i, pd = j, and {pg,pg+i) E E. In this case, the length of the path is d — 1. 
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An undirected graph is connected if for every pair of distinct nodes, there is a path connecting 
them. The Laplacian matrix of Q is defined hy L := A — A, where A is the diagonal matrix 
whose i-ih. diagonal is := ^ *^he symmetry of A). The incidence 

matrix B = [big] G of Q, with E := \£\, is that big := and bjg := for which 

(i, j) is the g-th edge of the graph (^f = 1,..., E), and big := 0 otherwise. By its construction, 
L = and B^ljsf = 0, where T and Ijv e R^ denote the matrix transpose and the vector 
of ones, respectively. See, e.g., IITSl for the details of graph theory. 

The stacking of vectors xi,... ,xn is denoted by [xi; • • • ; xn]- diag(i?i,..., Rn) is the block 
diagonal matrix with its Tth diagonal block Ri. The Kronecker product is denoted by $$ and B~^ 
denotes the Moore-Penrose pseudoinverse [|T^ of a matrix B. Euclidean norm is denoted by || ■ ||, 
while the 1-norm is denoted by || ■ ||i. For a compact set X, we define \\x\\x '■= minyg;^ ||x —j/||. 
Jjv is the identity matrix of dimension N and IT := - (l/iV)UlT. R>o denotes the set of 

nonnegative real numbers. 


II. Problem Formulation 

Fet us consider a group of N nonlinear systems, each of which is described by 

Xi f (^Xi^ Ui^ di]^ 

z = (1) 

Hi = h{xi), 

where x, G R" and Ui,yi G R'^ are the state, input, and output of the i-th system, respectively, and 
is incrementally output-feedback passive (iOFP) in the sense that there exist a storage function 
$ : R” X R"^ —)■ R>o, a number a G R, and two functions a and a of class /Coo such that 
a{\\xi - x'll) < $(xi,x9 < a{\\xi - x'|l) and 

x'i}f{xi, Ui, di) + x')/(x', m', <) (2) 

< - y[)^{yi - y'i) + (yi - y'i)^ ((«* + di) - («' + <)) 

hold for all {xi,Ui,di) G R” x R”? x R'^ and (x', «',(/') G R"" x R^ x R'^. Moreover, /(•,•,•) 
and h{-) are assumed to be locally Fipschitz and continuously differentiable in their arguments, 
respectively. The system satisfying (O is often referred to as cr-relaxed cocoercive system [@1. 
On the other hand, the system ([T]) is said to be incrementally passive if a < 0, and incrementally 
output-strictly passive if a < 0. Examples of iOFP systems include output-feedback passive 
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linear systems, Goodwin oseillators (see, e.g., [|71 Example 1]), and Van der Pol oseillators (see, 
e.g., SeetionlVT). 

The signal di represents the disturbanee acting on the z-th system and is generated by an 
exogenous system 

Wi = Si{wi), di = RiWi, Wi e (3) 

satisfying that Si(0) = 0 and for all {wi,w'j) G x 

{wi - {si{wi) - Si{w'i)) < 0. (4) 

Note that in this case, the solution of dH) always exists and any ball centered at the origin is 
forward invariant for ([3]). An example of such systems having the property (|4l) is the one given 
by Siiwi) = SiWi with Si skew-symmetric. 

We assume that the N systems dU) are defined over a connected undirected graph Q = 
with the node set M = {1,2,... ,N}. Each node z G AA is associated with the z- 
th system in dD and the edges in E define the interconnection structure. Under this setting, the 
objective is to design distributed controllers that enforce the asymptotic output synchronization 
of dB, i.e., 

lim \\yi{t) - yj{t)\\ = 0, i,j G J\f. 

t—>-oo 

The difference of our contribution from the vast majority of the literature on synchronization 
is that synchronization should be guaranteed despite the presence of exogenous time-varying 
signals which are different from node to node and which introduce a source of heterogeneity in 
the network. 


III. Controller Design for Output Synchronization 

In this note, we deal with two kinds of distributed controllers; one for controllers placed at 
the nodes (refer to, e.g., m-m or Section IIII-Ah with their inputs being 

p. = ^ay(z/j-z/i), z = l,...,iV, (5) 

and the other for controllers placed at the edges (see, e.g., [Utl, dUl or Section ITlI-BI) with their 
inputs constrained to be 

Sg — ^jgVj 1 


p = 1,...,E. 


( 6 ) 
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In both cases, the algebraie eonneetivity (say A 2 , the seeond smallest eigenvalue of L) of the 
graph Q, and internal model prineiple ifT^ play erueial roles for the synehronizability of the 
systems O as shown in BUl. Based on this observation, we provide internal-model-based 
solutions to the problem in the next subseetions, where the strong eoupling eondition A 2 > cr is 
no longer required. This relaxation is done by using one-dimensional adaptive eontroller at eaeh 
node tSeetion Illl-Ah or at eaeh edge tSeetion rill-BI) . 

A. Controllers at the nodes 

Homogeneous inerementally output-feedbaek systems of the form x* = /(x*) -t-Uj are known 
to synehronize for Ui = XljeA/'under suitable eonditions IIH. This is not guaranteed 
any longer if disturbanees di are affeeting the dynamies at the nodes. The eontrollers we propose 
in this seetion are a natural dynamieal extension of the statie diffusive eoupling, in whieh the 
dynamies is introdueed to deal with the heterogeneity of the disturbanees. 

Let us define x := und y ■= ^ Then, we have the following. 

Theorem 1. The outputs of the N systems o in closed-loop with the control 


Pi, 


(7a) 

ki 'liPi Pii 

li > 0, 

(7b) 

Ui Rifi T kiPi, 

i e M 

(7e) 


synchronize asymptotically whenever the closed-loop solution is bounded. In particular, lim^^oo WViit) — 
y{t)\\ = 0 holds for i E Af. 

Remark 1. A different approach to cope with the strong coupling condition is to use adaptive 
laws corresponding to the edges as in e.g., by employing, instead of ([5]) and (iTbl) . 

Pi = 'EjeN^ij^ijiyj ~ Pi) ^ij = Aijivi - yjYiyi - yj) with 

kij{0) = kij{0), yij = yji > 0 . ( 8 ) 

On the contrary, the update law (ITbl) is assigned to each node without the symmetry condition dS]). 

Thus, the number of adaptive laws in the entire system is less than that of 01/ , in general, 
yielding lesser dimension of controllers. In addition, only relative outputs of aggregated 
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into the variable pi in ([5]), are required to be measurable to implement (|7]), while additional 
communication of partial states of th internal models (TTal) is necessary in the case ofm. 

Proof: Let us consider the function Vi given by 

^i(^) = ^ xj), (9) 

ieAf jeAf 

where x := [a^i; • • • ',xn] £ We note that by Lemma [T] in the Appendix, there are p and 
rj of class /Coo such that r 7 (||x||) < Vi(x) < r 7 (||x||), in which x := (n 0 4) X = X — (Iat 0 x). 
Taking its time derivative along the solution of ([B and using Lemma [2] in the Appendix, we 
have 

Vi{x) < I X] X] aijiVi - VjViVi + - 4 + di- dj) 

i&M j&M ieAf jeAf 

= ay'^{L 0 Iq)y + y'^{L 0 Iq)u + y'^{L 0 Iq)d, 

where y := [yp • • • ; yjf\, u := [up • • • ; un], and d = [dp • • • ; d^]. Since u = —{KL®Iq)y — Ef, 

it further becomes 

4(x) < y~^{{aL - LKL) 0 Iq]y - y~^{L 0 Iq){Ri - d), 

where R := diag(4i,..., Rn), ^ := [4; • • • ; Civ], and iC := diag(fci,..., /cjv). 

Next, let us define := ^i — Wi and 4 := fcj — fc* with fc* G M to be determined, and consider 
the function 

ieAf i&Af ’ 

where ^ := [4;''' 'i^n] and k := [kp - ■ ■ ; fcAr]. Then, from the property dH), its time derivative 
becomes 

k) = Y^ (site) - SiK)) - 17 Rjpi + hpjpi} 

ieAf 

< -(T' p + p~^{k ® Iq)p = cF(L 0 Iq)y + y~^{LkL 0 Iq)y, 

where p := [pi; • • • ; pat], k := diag(fci,..., Izn) = K — k^N, and d := R^. 

We finally consider the Lyapunov function V{x, k) := Vi{x) + ¥ 2 ^, k), whose time deriva¬ 
tive is given by 

V(x, 7, k) < -y"' [{k^k - aL) 0 y. 


Let U = Q] be an orthogonal matrix sueh that U~^LU = A, where A = diag(Ai,..., Atv), 

Ai = 0, and Aj > 0 for i 7 ^ 1. Let us assume, with no loss of generality, that A 2 is the 
smallest nonzero eigenvalue of the Laplacian, and ehoose /c* > 0 suffieiently large to satisfy 
e := {KX 2 — a)X 2 > 0. Then, noting that QQ^ = If and If = 11^, one ean get 

V (x, i, k) < -ey'^iQQ'^ 0 Iq)y = -efy, (10) 

where yi := yi — y and y := [^ 1 ; ■ • • ;yAr] = (If 0 Iq)y. Integrating both sides of (fT^ . we 
have y{t) G £ 2 , i-e-, square integrable. Sinee the solution is bounded, Xi{t) and hence, yi{t) are 
bounded. The application of Barbalat lemma ifTSl . If20l proves the result. ■ 

In the proof of Theorem [H the function Vi{x) in ® is used to construct a valid Lyapunov 
function for the closed-loop system. If the system ([I]) admits a quadratic storage function, i.e., 
^{xi,x'i) = {xi — x'iyP{xi — x[) with P = P^ > 0, then by using Lemma |2] in the Appendix, 
we see that Vi{x) = x~^{L 0 P)x = 0 P)x which coincides with the one in ETl Theorem 

4]. Therefore, in the absence of d (because ll2T]l does not consider disturbances), our Lyapunov 
function can be considered as a nonlinear version of the one in ll 2 Tl Theorem 4]. 

As one can see from the proof of Theorem [U the design (|7]) provides some flexibility. When 
there is no external disturbance (i.e., ([U) is replaced by Xi = f{xi,Ui), yi = h{xi) and ^ holds 
with di = 0 and d\ = 0 ), the controller (TTbl) with its output Ui = kiPi, instead of dV]), solves 
the problem without requiring the strong coupling condition A 2 > a imposed in [|2|-[|4l. On the 
other hand, if A 2 > cr, then the use of the adaptive controller (iTbl) can be avoided. In this case, 
the controller is simply given by (TTal) and Ui = —Riii + Pi. 

B. Controllers at the edges 

In this subsection, we propose distributed controllers for synchronization, each of which is 
placed at the corresponding edge [HI, [|9]|. The input to the controller placed at the g-ih. edge is 
given by (jb]) and the input to the Ath system ([T]) is constrained to be 

E 

Ui ^ ^ ) i !)•••) R'y ( 11 ) 

9=1 

where Vg G is the output of the ^f-th controller. With this structural constraint, let us define 
H := (i?+ 0 Iq)R and let Hg G be the g-th block row of H = [Hi; ■ ■ ■ ;He], where 

R = diag(i?i, ... ,Rn) and m := This specific value of H in fact descends from the 


fulfillment of the regulator equations and the internal model property, shown in [[8|, whieh are 
neeessary to solve the output agreement problem for inerementally passive nonlinear systems 
under the eonstraint (fTTl) . In partieular, sinee u = —{B®Iq)v from (fTTl) and {B®Iq)Hw = ( 11 ® 
Iq)d from Lemma[3]in the Appendix, v = Hw (at steady-state) implies Ui+di = ^ 

(see (lH)), where v := [vi, ■ ■ ■ ;ve] and w := [wi; • • • ; wat]. In other words, the net effeet of 
disturbanees on eaeh system ([T]) is the same for all z = 1, ..., if v = Hw, whieh motivates 
us to eonsider eontrollers at the edges, able to generate the feedforward signal Vg = HgW. 

Theorem 2. The outputs of ([I]) in closed-loop with the control 


^iCg) + HJ Qg, 


(12a) 

^gQg Qg) 

o' 

A 

(12b) 

^9^9 ^gQg) 

g = l,...,E 

(12e) 


synchronize asymptotically if the corresponding closed-loop solution is bounded, where s(-) is 
defined by s{w) := [si(tei);--- ;sAr(teAr)]. Moreover, limf^oo ||2/j(f) — y(t)\\ = 0 holds for all 

Proof: Let us eonsider the funetion 

(13) 

i&N j&N 

whieh ean be shown to satisfy Z7(||5;||) < Wfix) < ^(||a;||) for some elass /Coo funetions rj and 
fj, again by Lemma [T| in the Appendix. Then, by using Lemma [2] with Oij replaeed by 1/N, one 
obtains the time derivative of Wi (x) as 

Wi{x) < ay^{If ® Iq)y + z/'^(n ® Iq)u + z/'^(n ® Iq)d. 

Let Q := [^i; • • • ; qe], C := [Cl ■ ■ ■ 'Xe\, H := diag(ifi,..., He), and K ;= diag(fi:i,..., ke). 
Noting that IfB = B, u = —{B ® Iq)v, and q = {B^ ® Iq)y, one further has 

Wfix) < ay'^{If ® Iq)y - z/^(SK5^ ® Iq)y + z/^(n ® Iq)d - z/^(S ® Iq)HC. 

Let us now define (g := (g — w and kg := Kg — where k* G M is a eonstant to be ehosen 
shortly. Let us eonsider the funetion 
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where ( := [Ci; • • • ; Ce] and k := [ki] • • • ; ke\- Then, using the property (01), we obtain its time 
derivative as 

E 

^2(C, «) = XI {C + H]Qg - s(w)) + KgolQg^ 

9=1 

< ~Ch'"q + 0 I,)q = 0 Ig)y + 0 I,)y, 

in whieh K ;= diag(Ki,..., ke) = K — 

Finally, let us consider the Lyapunov function W(x, k) = lLi(a;) + 1L2(C) for tho closed- 
loop system. Since ( = ( — {1e ® w) hy its construction and (11 0 Ig)d = {BB^ 0 Iq)Rw = 
{B 0 Iq)Hw = {B ® Iq)H (Ig 0 w) by Lemma [3] in the Appendix, the derivative of W{x, C, k) 
along the solutions of the system becomes 

W{x, C, k) < ay"'{Yi 0 Iq)y - n^y"'{L 0 Iq)y < - cr)y"^^, 

where A 2 and y are the ones in the proof of Theorem [T] If k* > 0 is chosen sufficiently large 
so that e := «:*A 2 — a > 0, then the result follows from the same arguments in the proof of 
Theorem [U ■ 

We note that the design (fT^ also possesses the flexibility discussed at the end of Section 
IIII-Al In particular, (I12bl) with Vg = KgQg solves the problem in the absence of disturbances and 
does not require the strong coupling condition in dUl, while (I12al) with Vg = Hg(^g + Qg achieves 
the asymptotic output synchronization of dB if A 2 > a. 

C. More on the proposed scheme 

Our scheme in the previous subsections can be further applied to the case in which the 
dynamics at the edges are given and not free to design ifTOll . ifTTll such as electrical networks [[22ll . 
[l23]l . An example of those cases is illustrated in Fig. [Tal where effective power transfer from the 
two sources to the load are aimed at. Such power transfer will occur when the terminal voltages 
of the two sources are synchronized or, equivalently, when the two nodes coupled through the 
dynamic edge in Fig. [Tb] are synchronized. Motivated by this example, let us consider a set of 
input strictly passive nonlinear systems, attached to the edges, 

% = '^gi.V9.Qg). Vg = VgijIg^Qg), C/ = 1, . . . , E, 








II 



Dynamic edge 



(a) An electrical network consisting of two voltage-dependent (b) An equivalent network of two nodes connected by a dynamic 


current sources and one load. 


edge, obtained from the Kron reduction process. 


Fig. 1: A motivating example of edges modeled by dynamieal systems. Y denotes the admittanee 
of the eleetrieal eireuit in the dotted box and vl is the voltage aeross the load. 


satisfying the dissipation inequality '^g{r]g) < 0g some positive definite '^g, where 

%ljg and ifg are loeally Lipschitz and eontinuous, respectively. 

We now provide the following result. 

Proposition 1. Assume © is disturbance-free, i.e., Xi = f{xi,Ui), yi = h{xi), and satisfies Q 
with di = d[ = 0. Then, their outputs synchronize asymptotically if Ui is given by (fTTl) . the 
corresponding closed-loop solution is bounded, and the strong coupling condition A 2 > cr is 
satisfied. 

Proof: Let us consider the storage function W{x, rj) = Wi{x) + X]f=i ^^giVg)^ where Wi{x) 
is in (fT3l) and r] := [rji]--- ;//£:]. Then, noting u = —{B 0 Iq)v and g = {B~^ ® Iq)y, its time 
derivative is given by 

W{x, 7]) < ay^{Il 0 Iq)y + ?/"^(n (g) Iq)u - g^g + g 

= ay"'{Yi (g Iq)y - y"'{BB"' g Iq)y < -(A2 - 

The rest can be proven by following the similar arguments in the proofs of Theorems [U and [2l 

■ 

The result of the proposition finds its application in problems of synchronization of, e.g., 
(a) electrical networks without shunt elements and with dynamic heterogeneous edges and (b) 
electrical networks with shunt elements in which their Kron-reduced networkil] contain identical 


'Refer to, e.g., (24l for the details on the Kron reduction process of electrical networks. 
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shunt elements with (possibly) heterogeneous dynamie edges and eaeh node eonsisting of a 
souree and its eorresponding shunt element in the redueed networks is iOFP (refer to Fig. [T] 
as an example of this ease). Therefore, Proposition [T] offers a eomplementary result and proof 
teehnique to If22l . [l23ll . We note however that in all eases, boundedness of the elosed-loop 
solutions is indispensable for proving the synehronization. Sueh boundedness and a elass of 
systems ineluding, e.g.. Van der Pol oseillators will be diseussed in Seetion |IVl Oseillators 
whieh are not globally Lipsehitz sueh as Van der Pol oseillators were not handled in [l22ll . [|2^ . 

When there are disturbanees aeting on the systems, the problem beeomes more ehallenging. 
However, at least on a eomplete graph, the outputs of the systems ean be synehronized as follows. 

Proposition 2. Suppose that Oij = a for i,j E M and for some a > 0. Then, the outputs of the 
N systems ([I]) in closed-loop with the controllers 

E 

Pii tli ^ ^ ^ig^gy i — 1) 2, • • • , N (14) 

9=1 

synchronize asymptotically if X 2 > cr and the corresponding closed-loop solution is bounded. 

Proof: Let us eonsider the storage funetion W{x, p, f) = ^7ii/ 

{2aN), where := ii — Wi and ^ ;= [^ 1 ; • • • ; ^n]. Taking its time derivative and bearing in mind 
that p = —aN{Il 0 Ig)y and u = —RS, — {B ® Ig)v, we get 

w{x, p,7) < a?/"^(n (g) Iq)y + |/'^(n ® Iq){u + d) - g + v~^g - (Fp/{aN) 

= -y^ (L - o-n (g) Ig) y < -(A2 - a)y~^y, 

where d := Rf — d. Thus, the result again follows from the similar argument in the previous 
subseetions. ■ 

It is also possible to aehieve the output synehronization without having the strong eoupling 
eonditions imposed in the previous propositions. The idea is to eonsider the passive dynamieal 
edges 

Vg = ^g{VgyQg), Vig = ipg{pg, Pg), g = l,...,E, 

satisfying ^g{pg) < vjgpg and to assign adaptive laws to the eorresponding edges. We note that 
sueh assignment may not be feasible in the ease of, e.g., eleetrieal networks. However, the two 
eorollaries given below are theoretieally interesting and eomplement the results of the previous 
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propositions in the sense that they provide adaptive variants of the propositions to relax strong 
eoupling eonditions. 

Corollary 1. Assume the setup of Proposition [7] Then, the outputs of ([U) synchronize asymp¬ 
totically by the control 

kg ^gQg Qgi ^2g ^gQgy 9 ■ ■ ■ y ^ y 

E 

Mi = - ^ hig {vig + V2g) , i = 

ff=l 

if the corresponding closed-loop solution is bounded. 

Proof: Let us eonsider the funetion W (x, p, k) = Wi{x) + ^gidg) + X]f=i 

where kg := Kg — > 0, and k := [kp, • • • ; Mg]. Let us define m* := [vn] ■ ■ ■ ',ViE] with 

i = 1, 2, and v = vi V 2 . Then, noting u = —{B ® Iq)v and M 2 = (K (g) Iq)g = {KB^ (g) Iq)y 
with K := diag(Ki ,... ,ke), one obtains the time derivative of W as 

W{x,p,k) < ay^{U® Iq)y{U® Iq)u-\-vjg-k g^{k® Iq)g 
= ay'^iU (g) Iq)y -y'^{B® Iq)v2 + y"'{BKB^ ® Iq)y 
= ay^{Yi ® Iq)y - K^y^{BB^ ® Iq)y < -(k*A2 - (T)y~^y. 

The rest ean be proven by following the similar arguments in the proofs of Theorems [U and [2l 


Corollary 2. Assume the setup of Proposition \2i Then, the outputs of the systems ([B in closed- 
loop with 


Pi, kg ^gQg Qgy ^ 0, 


E 

Ui Ri^i ^ ^ big iv\g T '^2g') y ^2g 
9=1 


synchronize asymptotically whenever the solution is bounded. 


Kggg 


Proof: Let us eonsider the funetion W{x, p, k) = Wi{x)-kYl^=i ^ g{'9g)~kYl!i=i C^6/(2a-/V) + 
X]f=i «g/(25g)- Then, noting u = -Rf-{B®Iq)v, p = -aN(JI®Iq)y, and V 2 = {KB~^ ®Iq)y, 
we have 

W{x, p, k) < ay^ijl (g) Iq)y + p'^(n (g) Iq){u + d) + vjg - ^p/{aN) + p'^(K (g) Iq)g 
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= (Ty'{U ® Iq)y - y^{B 0 Iq)v2 + y^ {BKB^ ® Iq)y < -(k*A2 - (T)y'^y. 

Thus, the result again follows from the similar arguments in the previous subsections. ■ 

The extension of these results to the case of dynamics at the edges satisfying different 
dissipation inequalities from those assumed previously and to the case of graphs which are 
not complete, as well as the relaxation of the strong coupling condition by assigning adaptive 
laws to the nodes, is left for future investigation. 

IV. Solution Boundedness 

In the previous section, the synchronization is guaranteed, provided that the solution of the 
closed-loop system is bounded. We discuss in this section on what conditions of the open-loop 
system ([T]) such boundedness is ensured under the control (|7]) or (fT^ . To do this, we further 
assume that the output map h{-) is globally Lipschitz and there is a compact set X C M”, 
invariant for ([1]) when = 0 and di = 0, such that the system ([U) is input-to-state stable (ISS) 
with respect to Ui := [ui]di] relative to X. In other words, there are two functions j3{-, •) and 
y{-) of class ICC and /C, respectively, such that the solution of ([B satisfies 

/^(Iki(O)IU,^),/^ f sup ||Mi(r)||^l. (15) 

\0<T<t J ) 

We refer the reader to [l25l Section III] for some details of the condition (fT5l) . A class of systems 
ensuring both conditions Q and (fTSl) will be discussed after presenting the following result. 

Proposition 3. Let the ISS property (fTSl) be satisfied and suppose either ([7]) or (fT2l) is ap¬ 
plied to the system o. Then, the solutions of the closed-loop system are bounded and satisfy 
limt^oo Wdiit) - y{t)\\ =0, i = 

Proof: We prove the case for controllers placed at the nodes (|7|)- Other cases in Sections 
IIII-BI and IIII-CI can be proven similarly. 

Let the variables x and C, be the ones in the proof of Theorem [T1 and define k := [ki] - ■ ■ ^k^f. 
Let [0,T,i), Tu < -l-C )0 be the maximal time-interval, where the unique solution of the closed- 
loop system starting at (x(0), .^(0), fc(0)) exists. Then, (fTOl) holds on this interval. As a result, 
^{t) and k{t) are bounded for all t G [0,T„) because w{t) is bounded for f > 0 from dH). 
Note that xfit) — Xj{t), i,jeJ\fare also bounded on the interval since Vi{x) in dH) is positive 
definite with respect to x. Together with the globally Lipschitz output map h{xi), this implies 


\xi{t)\\x < max 
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Fig. 2: A Lur’e-type nonlinear system eonsisting of a passive linear system E and a statie 
nonlinearity 0(-) in the feedbaek path. 


the boundedness of the eontrol inputs Ui{t), i G TV on [0,T„). Therefore, the solution Xi{t) is 
bounded from the ISS property (fTSl) . leading to the existenee of a positive eonstant M sueh 
that \\[x{t)-, ^(t); k{t)]\\ < M for all t G [0,T„). This means that the solution of the elosed-loop 
system ean be extended to [0,T] for some T > T^, whieh eontradiets to the assumption that 
[0,T,i) be the maximal interval of the solution existenee. Thus, the result follows if we repeat 
the previous proeess with T„ = +cx). ■ 

Theorem [T] (or ^ eombined with Proposition [3] requires the open-loop system ([T]) to satisfy 
Q and (fT5l) . simultaneously. We note that sueh systems indeed exist. To see this, let us eonsider 
the Lur’e-type nonlinear systems (see Fig. [2] or ifTSl Seetion 7]) of the form: 

Xi = Axi + B + Ui + di), yi = Cxi, (16) 

where (C, A) is deteetable and the linear part is passive, i.e., there is a matrix P = > 0 sueh 

that A~^P -f PA < 0 and PB = C~^ hold. The nonlinearity 0 : M — >■ M is loeally Lipsehitz and 
satisfies that lim.r^+oo = +oo and limT-^_oo 0(r) = —cxd. It is also assumed that there is 
r* > 0 sueh that 0(r) is monotonieally inereasing on both intervals (— cxo, —r*] and [r*, -f cx)J]. 
Then, we have the following. 

Proposition 4. The Lur’e feedback system (fT^ satisfies both of the iOFP condition dH) and the 
ISS property (fT5l) . 

Proof: The satisfaetion of the ISS property (fTSl) follows from Il27l Thoerem 2]. To show 


^An example of such nonlinearities is the function <}>{■) illustrated in Fig.|2 We note that in contrast to 1261 . the nonlinearity cj> 
considered here is neither incrementally passive nor incrementally sector hounded due to (possible) negative slope at the origin. 
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the iOFP of (fT^ . let us consider a storage function = ^{xi — x'j)~^P{xi — x[) and its 

time derivative along the solutions Xi and x' of (dH) as 

^{Xi,x[) = {xi - x'iyP{Axi + B{-(j){yi) + u* + di) - Ax[ - + u[ + <)} 

< iVi - y'i) (- 0 ( 2/0 + 0 ( 2 /D +Ui-u'i + di- d'i). 

This implies that (fT^ is iOFP if there is a > 0 satisfying 7r(j/j, y'^) := —{yi—y'i) { 4 >{yi) — 4 >{y'd) ^ 
a{yi — y'iY- The existence of such a can be shown as follows. 

i) yi,y'i e [—r*,r*]: Since 0 is locally Lipschitz, vr < w{yi — y[Y holds, where cu > 0 is a 
Lipschitz constant of 0 on [—r*,r*]. 

ii) yi, y[ G [r*, +cx)): By the monotonicity of 0 on [r*, +oo), it holds that vr < 0 < w{yi — y[Y. 

iii) yt G [—r*,r*] and y[ G [r*,+oo): From the Lipschitz property and monotonicity of 0, it 
follows that 

TT = -{yi - y'i) (0(2/0 - 0 (t-*)) - (2/i - y'i) ( 0 (n) - 0(2/')) 

< ^\yi - y'iWyi - < ^{yi - y'if- 

iv) yi G (— 00 , —r*] and ?/■ G [r*, +cxo): Again, by the Lipschitzness and monotonicity of cj), we 
have 

TT = -{yi - y'i) (0(2/0 - 0(-a)) - {yi - 2 /O (0(-a) - 0(r*)) - {yi - y'A {(j){T^) - 0(i/')) 

< ^\yi - 2/(11 (-a) - r*| < w{yi - y^f. 

The remaining cases can be proven similarly. Therefore, TT{yi,y'i) < a{yi — y'iY holds for all 
yi, 2 /( G M with a = w > 0. ■ 

Examples of such Lur’e-type nonlinear systems include Van der Pol oscillators (see Section 
|V]for the details) and the Lienard-type dead-zone oscillators (refer to ll23ll ). Another example is 
the Chua’s circuit If28l , Il29ll whose dimensionless form is given by 

ii,l = Ci{Xi^2 — Xi^i — (j){Xi^i) + Ui + di), 

Xi,2 = Xi^i - Xi^2 + a;*,3, (17) 

(^2,3 Vi 

where x* := [xiX-, Xi^ 2 ', Xi^A ^ ci,C 2 > 0, and the nonlinearity 0(0 is shown in Fig. [3l It is 
noted that the unforced system (i.e., (fTTl) with u, = 0 and di = 0) with particular selection of 
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Fig. 3: The characteristic of the nonlinearity ((){■) in Chua’s circuit. 


parameters Cj, n, r*, and zj has the double-scroll (chaotic) attractor (see, for instance, [l28l Fig. 
6] and [|29l Fig. 23]). The Chua’s circuit (fT71) is of the form (fT^ and its linear part satisfies the 
passivity requirements P + PA < 0 and PB = C'^ with P = diag(l/ci, 1, 1 /C 2 ). Moreover, 
the piecewise linear function 0(-) in Fig. [3] guarantees all the required properties such as local 
Lipschitzness, monotonicity, and limiting behaviors at the infinity. Therefore, by Proposition 01 
the Chua’s circuit (fTTl) ensures both of the iOFP and ISS conditions. 

V. Computer Simulation: Van der Pol Oscillators 
Let us consider Van der Pol oscillators [[l5l Example 2.9] given by 

^i,2i 

Xi^2 =xl2/‘i)+Ui + di, 7 = 1,..., 4, (18) 

Vi = Xi^2, 

where Xi := [xi^i;Xi^ 2 ] ^ and i/ G M is a positive constant. If we define A := [0 1; —1 0], 
B := [0; 1], C := [0 1], and (j){yi) := z/(a;f 2/3 —then (fTSl) is of the form (fT^ and satisfies all 
the required properties with P = I 2 . Thus, the Van der Pol oscillator (fTSl) is iOFP and satisfies the 
ISS property given in (fTSl) . Moreover, one can verify that a > u. We set u = 1. Meanwhile, the 
disturbances dj, 7 = 1,..., 4 in (fTSl) are assumed to be generated by the exosystems (satisfying 
the requirement (0])) 

di TUz, 


Wi = 0, 


TUj G M, * = 1, 2, 
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(a) Trajectories of the states 


(b) Trajectories of the states 



(c) Trajectories of the outputs yi{t) = Xi^ 2 {t). 



(d) Trajectories of the outputs yi{t) = Xi^ 2 {t). 


Fig. 4: Simulation results with controllers placed at the nodes (left) and at the edges (right). 


Wi = Wi, di = 



1 0 

RiW-i 


Wi, Wi G i = 3,4. 


(19) 


In this setting, we first perform a set of simulations to demonstrate the results in Sections 
IIII-AI and ITTT-BI with the interconnection structure Q characterized by its Laplacian L = 0.09L 
and incidence matrix B = 0.35, where 



' 5 

-1 

0 

- 4 ' 
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0 

0 

0 
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-1 
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-9 

10 

-1 
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-3 

0 

1 

0 


-4 

-4 

-1 

9 


0 

0 

-2 

-1 

-2 


Note that the strong coupling condition is not satisfied in this case since A 2 = 0.4 and cr > 1. 

































19 


Figs. andshow a simulation result with the eontrollers plaeed at the nodes (|7]), while Figs. 
|4b] and Hd] eorrespond to a simulation result with the eontrollers at the edges (fT^ . The update 
gains 7 j and 5g are seleeted as 7 * = 5^ = 1, z = 1,..., 4, = 1,..., 5. All the elements of the 

initial eonditions of the Van der Pol oscillators and exosystems are randomly chosen within the 
interval [—3, 3], while those of controllers are set to all zeros. From Figs. S^and |4dl it is seen 
that the outputs of the Van der Pol oscillators (fTSl) synchronize asymptotically as expected. Note 
however that the remaining states z = 1,... ,4 do not synchronize since the proposed 

controllers guarantee the output synchronization only. 

For the case of edges with dynamic systems attached as in Section ITTT-Cl let us consider the 
systems (flSl) connected over the complete graph with = 1 for z, j = 1,..., 4, i.e., 




' 1 

1 

1 

0 

0 

0 



-1 

0 

0 

1 

1 

0 

L = 4n, 

B = 

0 

-1 

0 

0 

-1 

1 



_ 0 

0 

-1 

-1 

0 

-1 


In this case, the strong coupling condition is guaranteed since A 2 = 4 and the shortage of 
passivity, a, can be taken as a = iz = 1. The edges connected to the first Van der Pol oscillator 
are modeled by 

% = ^ 5 , Vg = 'ng + Qg, 5(=1,2,3, 

whereas the rest are given by the dynamical systems 

Vg = -Vg + Qg, Vg = T]g + Qg g = A, 5, 6. 

Note that with the storage function '^{rig) = \rfg, both dynamics satisfy the dissipation inequality 
for the input strict passivity given in Section ITTT-C[ The left column of Fig. [5] shows a simulation 
result with the controller (fTTI) in the absence of disturbances, while the right column of Fig. [ 5 ] 
shows a simulation result with the controller (fT4l) in the presence of the external disturbances 
generated by the systems (fT9l) . Each element of the initial conditions of the Van der Pol 
oscillators, exosystems, and dynamics at the edges are randomly chosen within the interval 
[—3,3], while those of controllers are set to all zeros. As expected from Propositions [H and [2l 
the outputs of the Van der Pol oscillators (fTSl) synchronize asymptotically. See the middle row 
of Fig. [51 i.e.. Figs. [5c] and [5dl 
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(a) Trajectories of the states 


(b) Trajectories of the states 




(c) Trajectories of the outputs yi{t) = Xi^ 2 {t). 


(d) Trajectories of the outputs yi{t) = Xi^ 2 {t). 




(e) Trajectories of the edge states r]g{t). 


(f) Trajectories of the edge states 'r]g(t). 


Fig. 5: Simulation results in the presenee of dynamieal edges without disturbanees (left) and 
with disturbances (right). 
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VI. Conclusion 

In this note, we have eonsidered the asymptotie output synehronization problem of inere- 
mentally output-feedbaek passive nonlinear systems (defined over a eonneeted undireeted graph) 
in the presenee of external disturbanees. Two different struetures of distributed eontrollers are 
eonsidered; one plaeed at the nodes and the other plaeed at the edges. In both oases, the proposed 
eontrollers are shown to solve the problem if the solution of the eorresponding elosed-loop system 
is bounded. A pair of results has been also disoussed to deal with the ease in whieh the dynamios 
at the edges are given. A olass of inorementally output-feedbaek passive systems that guarantee 
the boundedness of the elosed-loop solutions is then provided, by using the notion of input- 
to-state stability relative to a set. The eontrollers are synthesized based on the adaptive eontrol 
teehnique and the internal model prineiple and, as a eonsequenee, do not require the so-ealled 
strong eoupling eondition. 


Appendix 

Lemma 1. Suppose that Q is connected. Then, there are two class /Coo functions rj and fj that 
satisfy ? 7 (||a;||) < Vi(a;) < ?/(||5i||), where Vi{x) is in dH) and x = (n® 

Proof: (Existence ofr]): Define oC'i •— ^2 oC • Sinee 

a(ri H-h tn) < a{NTi) H-h a{NTN) (20) 

holdj^ for any elass /Coo funetion a, one has 

Vi{x) < + 

= ^Z\ia(2||xi||) < Z\max^a(2||xi||), 
ieV i&M 

where A^ax := max* > 0. Defining X := ; H^atII] G the inequality further 

beeomes 

Vi{x) < a(2||xi|| + • • ■ + 2||x^||) = NA^^a{2\\X\\,) < NA^,Jy{2^/N\\X\\). 

ieJV 

^Indeed, letting Ti be one of the largest elements among ri,..., tjv , we have that Q:(ri -|- • • • -I- tn) < a{NTi) < a{NTi) + 
■ • • -I- a{NTN). 
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Observing ||X|| = ||x||, set r 7 (r) := iVZ\max«(2\/iVr). 

(Existence ofr]): For I <i < N — 1 and j > i, let {pij,i,Pij, 2 , ■ ■ ■ iPijAj} be one of the shortest 
paths eonneeting the nodes i and j. Its length is dij — 1. The total length of those N{N — l)/2 
paths satisfies 

N{N-lf 


N-l 


N -1 




2=1 j>i 2=1 j>i 

Define Omin := mina-^.^o > 0 and let A/i be the set of neighbors of node i. Then, using 
and triangular inequality, and noting that Pij^i = i and Pij^dij = j, we have 

TV—1 dij — 1 


Vi{x) > Y - ^ill) ^ 


i&N j&Ni 


^min 

NiN -If ^ 

^ 2=1 j >2 g=l 




TV—1 / — l II 

Omin ^ ^ ^ j 2-^g=l \\^Pij,g ^Pij, 


> 


9 + 1 I 


N(N -ly 

^ ' i=\ jyi 


dij 1 


N-l 


> 


EEa 


I X i X j 

N-l 


N N 




X i X i 


N(N -ly ^ N-l J 2N(N -ly ^ N-l 

'' > 2=1 j>i V / ^ ^ 2=1 j = l V 

By using (l20l) and triangular inequality again, one finally has 


TV TV 


vyx) > 


2N(N - 1)2- yN(N - 1)2 ^ ^ 


OC 7 OC 7 


N 


> 


rtt 


3712^11**11 T: 


2Ar(iV- 1)2- UiV- 1)2 ^ 
Set p{t) := a ^^M r/(N - iy)/{2N{N - 1)2}. 


a 


ll^ll 


2Ar(iV- 1)2- \{N-1) 


Lemma 2. Let Oi^Di G for i = 1,..., A^. Then, the following holds. 

} E Z “**('’* - ® 4)'^. 

i&M j&M 

where 6 := [ 6 ^ 1 ; ■ ■ ■ ; 9n] and -d := [ t ^ i ; • • • ; 1}^]. 


Proof: The following eomputation proves the lemma. 

7 E E “«('** - = E f E (A ® /,)tf +0 

i&N j&M i&M i&M 


= 6' (A® Ig)'d - 6'^ (A ® Ig)d = e'(L® Ig)'d, 


T/ 


where Ai is the i-th row of symmetrie adjaeeney matrix A. 
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Lemma 3. Let Q be a connected undirected graph. Then, BB^ = LL^ = 11 holds. 

Proof: Since B^ = B^{BB^)^ by IIT^ p. 49] and L = BB^, we have BB^ = LL^. On 
the other hand, by ifT^ p. 60, Corollary 7], LL+ is the unique (orthogonal) projector on im(L) 
along ker(L+) = ker(L^) = ker(L), where im(L) and ker(L) denote the image and kernel of 
L, respectively. Therefore, LL+ = If because 11 is also the projector on = im(L) along 
span(lAr) = ker(L+), where is the orthogonal eomplement of Itv- ■ 
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